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ABSTRACT 

Context. Strong lensing is one of the most direct probes of the mass distribution in the inner regions of galaxy clusters. It can be used to 
constrain the density profiles and to measure the mass of the lenses. Moreover, the abundance of strong lensing events can be used to constrain 
the structure formation and the cosmological parameters through the so-called "arc-statistics" approach. However, several issues related to the 
usage of strong lensing clusters in cosmological applications are still controversial, leading to the suspect that several biases may affect this 
very peculiar class of objects. 

Aims. With this study we aim at better understanding the properties of galaxy clusters which can potentially act as strong lenses. 
Methods. We do so by investigating the properties of a large sample of galaxy clusters extracted from the N-body/hydrodynamical simulation 
MareNostrum Universe. We perform ray-tracing simulations with each of them identifying those objects which are capable to produce strong 
lensing effects. We explore the correlation between the cross section for lensing and many properties of clusters, like the mass, the three- 
dimensional and projected shapes, their concentrations, the X-ray luminosity and the dynamical activity. 

Results. We quantify the minimal cluster mass required for producing both multiple images and large distortions. While we do not measure a 
significant excess of triaxiality in strong lensing clusters, we find that the probability of strong alignments between the major axes of the lenses 
and the line of sight is a growing function of the lensing cross section. In projection, the strong lenses appear rounder within i^oo, but we find 
that their cores tend to be more elliptical as the lensing cross section increases. As a result of the orientation bias, we also find that the cluster 
concentrations estimated from the projected density profiles tend to be biased high. The X-ray luminosity of strong lensing clusters tend to be 
higher than that of normal lenses of similar mass and redshift. This is particular significant for the least massive lenses. Finally, we find that 
the strongest lenses generally exhibit an excess of kinetic energy within the virial radius, thus indicating that they are more dynamically active 
than usual clusters. 

Conclusions. We conclude that strong lensing clusters are a very peculiar class of objects, affected by many selection biases which need to be 
properly modeled when using them to study the inner structure of galaxy clusters or to constrain the cosmological parameters. 

Key words. Keywords should be given 



1. Introduction 

Gravitational lensing is one of the most powerful tools for 
studying the formation of cosmic structures in the universe. 
The light from distant sources, traveling in space and time, is 
deflected by the matter along its path before being collected 
by the observers. Thus, we measure an integrated effect which 
contains a wealth of information about the cosmic structures at 
different epochs. 

Depending on the impact parameter on the intervening mat- 
ter and on the mass of the deflectors encountered by the light 
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along its path, gravitational lensing manifests itself in the weak 
and in the strong regimes. In the weak lensing regime the 
shapes of distant galaxies, which happen to be at large angu- 
lar distances from the largest mass concentrations on the sky, 
are slightly changed, such that this effect can be revealed only 
though statistical measurements. Nevertheless, these tiny dis- 
tortions can be used for tracing the large scale structure of the 



universe both in two and in three dimensions (see e.g. Fu et al. 
[2008} |Benjamin et al.|2007| |Massey et alT1[20Tj7l for some re- 
cent results), from which important cosmological constraints 
can be derived (Bartelmann & Schneider 2001 ). This is a field 
of research which will have extraordinary improvements in the 
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next decades, thanks t o some upcoming missions (|Wittman| 
et al.|2006||Kaiser|2007||Jelinsky & SNAP Collaboration|2006| 
Refregier et al.||2008| ). Weak lensing allows to reconstruct the 



|2007al|Kneib et al.|2003irGavazzi et al.|2003"l|Broadhurst et al.| 



mass distribution up to the outskirts of galaxy clusters (see e.g. 
|Dahle|20061 [Clowe et al.|2006l|Hoekstra|2007l ). 

Strong lensing is an highly non-linear and relatively rare 
effect which is observable in the central regions of galaxies and 
clusters. In this regime, the background sources can be mul- 
tiply imaged and/or highly distorted to form very elongated 
images, the so called "gravitational arcs". They are powerful 
cosmological probes for many reasons. First, such events can 
be used to investigate the inner regions of the lenses. Thus, 
they can be used to test the predictions of the Cold-Dark- 
Matter paradigm on the inner structure of dark matter halos 
(|Meneghetti et al.|[200T) |Mao et al.||2004| |Sand et al.|[20( 
B artelma nn"& Meneghetti| |2004| |Meneghetti et al.||2007b 
Limousin et al.|2008 1. Second, they can be used to recover the 



mass distribution in the centre of the lenses, providing com- 
plementary informations to those obtained from weak lensing 



(Bradac et al.||2005| Diego et al. 


2005| Cacciato et al. 2006| 


Limousin et al. 2007 Merten et a 


. 2009 ). Lensing masses can 



then be used for measuring the cluster mass function. Third, the 
position and the distortions of the strongly lensed images as a 
function of the source redshift reflect the geometry of the uni- 
verse ( |Soucail et al.|2004| [Meneghetti et al.|2005| >. Finally, the 



probability of observing strong lensing events is deeply con 
nected to the abundance, the mass and the formation epoch 
(through the concentration) of the lenses. This makes statis- 
tical lensing a potentially powerful tool to study the structure 
formation ( |Bartelmann et al.|1998) |Li et al.|2005j [Meneghetti 



et al. 2005) 

In this paper, we focus on the properties of the most mas- 
sive and therefore most efficient strong lenses in the universe: 
the galaxy clusters. In the framework of the hierarchical sce- 
nario of structure formation, these are the youngest bound sys- 
tems in the sky. About 85% of their mass is believed to be 
in the form of cold-dark-matter (Gottlober & Yepes 2007 ). 
The remaining 15% is made of a diffuse gas component, the 
Intra-Cluster-Medium, and of other baryons in the form of 
stars, the vast majority of which is inside the cluster galaxies. 
Being relatively young structures, the interaction between the 
baryons and the dark matter is less strong than in older systems 
like galaxies. For this reason, clusters are important laborato- 
ries for studying the properties of the dark matter (Markevitch 
et al.|[2004] ). However, there are several issues that we need 
to take into account when studying these systems. In partic- 
ular, clusters where gravitational arcs are observed are a lim- 
ited fraction of the total number (|Sand et al. |2005 Luppino 
et al.||1999[|Gladders et al.||2003[ |ZanTsky & Gonzalez ||2003' 



Hennawi et al. 2008 1 and therefore a particular class of ob- 
jects. Broadly speaking they are the most massive clusters, but 
there are several other properties that boost the cluster abil- 
ity to produce strong lensing events. For example, we know 
that substructures, asymmetries and projected ellipticity of the 
lenses are all contributing to the strong lensing cross section 
of clusters ( |Meneghetti et al.|2003"b} |Meneghetti et al.|2007a) . 
Both observations and simulations agree that strong lenses 



2008 ). Triaxiality is also relevant, because clusters seen along 



their major axis are more efficient lenses (Ogu ri et al.||2003 



Although cluster galaxies statistically do not change the dis- 
tributions of the arc properties ( Meneghetti et al.|2000[ |Flores 



|et al.|2000| |Hilbert et al.|2008| ), cD galaxies sitting at the bot- 
tom of the cluster potential well increase the ability for pro- 
ducing long and thin arcs by ~ 30 - 50% (Meneghetti et al. 



2003a). The gas physics, in particular cooling, could also affect 



the strong lensing properties of clusters ( Puchwein et al.|2005 



Wambsg anss et al.|2008[ |Mead et al.|2010| ). Finally, |Torri et at 
(2004} showed that the cluster ability to produce gravitational 
arcs can also be enhanced by the dynamical activity in the lens. 
By studying with high time resolution how the lensing cross 
section changes during an edge-on collision between the main 
cluster clump and a substructure, these authors found that the 
strong lensing efficiency is boosted by a factor of 10 during the 
merging phase. Later, Fedeli et al. ( 2006|), using semi-analytic 



methods, showed that the arc optical depth produced by clus- 
ters with moderated and high redshifts is more than doubled by 
mergers. 

Although an extensive work has been done in the past 
decade, a better characterization of the strong lens cluster pop- 
ulation is mandatory. Given the complexity of clusters and 
the great importance that many of their properties have for 
strong lensing, the only reliable way to do that is through 
the ray-tracing analysis of a large number of simulated clus- 
ters. A first important work in this framework was done by 



Hennawi et al. (2007i, who analyzed a sample of 878 clus- 
ters from an N-Body pure dark-matter cosmological simula- 
tion. Important properties like concentrations, axis ratios and 
substructures were discussed. In this work, we aim at extend- 
ing the analysis of Hennawi et al. (2007]) in three ways. First, 



have high concentrations (Hennawi et al. 2007 Fedeli et al. 



we include a much larger number of clusters (now 49366 sys- 
tems), taken from a larger cosmological volume (500/i~'Mpc 
vs 320/i~'Mpc). Second, the clusters used here are obtained 
from an N-body-hydrodynamical simulation where the evolu- 
tion of the gas component is also considered. Thus, we can 
correlate the lensing properties of clusters with some important 
X-ray observables. Third, we study in detail the possible corre- 
lations between cluster dynamics and strong lensing efficiency, 
which was made only through analytical models so far. 

The plan of the paper is as follows. In Sect. 2 we sum- 
marize the main characteristics of the cosmological simulation 
MareNostrum Universe. In Sect. 3 we discuss the simulation 
methods and we define several lensing quantities useful for the 
following analysis. In Sect. 4 we discuss the correlation be- 
tween lens masses and strong lensing ability. Sect. 5 is dedi- 
cated to the statistical analysis of the shapes and orientations 
of strong lensing clusters. We discuss the biases in the concen- 
trations in Sect. 6. In Sect. 7 we focus on the X-ray proper- 
ties of strong lensing clusters. Finally, in Sect. 8 we correlate 
the strong lensing efficiency with the dynamical state of the 
lenses. We summarize the main results and the conclusions of 
this study in Sect. 9. 
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2. The MareNostrum Universe simulation 



The MareNostrum Universe (Gottlober & Yepes 2007) is a 
large scale cosmological non-radiative SPH simulation per- 
formed with the Gadget2 code (Springel 2005). We briefly 
summarize the relevant characteristics here and we refer the 



reader to the paper by Gottlober & Yepes (2007 ) for a more de- 
tailed description of the simulation. This was run in 2005, dur- 
ing the testing period of the MareNostrum supercomputer us- 
ing the WMAP1 normalisation, namely £2 m ,o = 0.3, Qa,o = 0.7 
and erg = 0.9 with a scale invariant primordial power spec- 
trum.. After the release of the 3-year WMAP data the simula- 
tion has been repeated at lower resolution with the predicted 
low normalization ( Spe rgel et al.|2007" i and with a higher nor- 
malization of cr 8 = 0.8 which is better in agreement with the 
5-year WMAP data ( |Komatsu et al.||2009) . Comparing pre- 
dictions from these numerical simulations with recent obser- 
vational estimates of the cluster X-ray temperature functions 



Yepes et al. (2007} argue that the low normalization cosmo- 



logical model inferred from the 3 year WMAP data results is 
barely compatible with the present epoch X-ray cluster abun- 
dances. Therefore, the original WMAP1 normalized version of 
the simulation is used in this paper. The simulation consists of 
a comoving box size of 500 hr l Mpc containing 1024 3 dark 
matter particles and 1024 3 gas particles. The mass of each 
dark matter particle equals 8.24 x 10 9 M e h~ l , and that of each 
gas particle, for which only adiabatic physics is implemented, 
is 1.45 x 10 9 M Q h~ l . The baryon density parameter is set to 
£2b,o = 0.045. The spatial force resolution is set to an equiv- 
alent Plummer gravitational softening of 15 h -1 kpc, and the 
SPH smoothing length was set to the 40th neighbour to each 
particle. 

In order to find all structures and substructures within the 
distribution of 2 billion particles and to determine their prop- 
erties we use a hierarchical friends -of-friends (FOF) algorithm 
( |Klypin et aL"1|1999| l. With a basic linking length set to 0.17 
times the mean interparticle distance, we extract the FOF ob- 
jects at all redshifts. The final catalog of identified objects con- 
tains more than 2 million objects with more than 20 DM par- 
ticles at z = 0. The same objects and their progenitors are 
contained in the catalogs corresponding to higher redshift out- 
puts of the simulation. In this sense a correlation exists between 
different redshift slices. In a second step we divide the linking 
length by 2" (n — 1, 3) to find substructures of the clusters. In 
particular we use n — 2 to identify the highest density peak 
which we associate with the center of the cluster. 

All the FOF groups with mass larger than 10 13 hr l M Q are 
then stored into sub-boxes of cubic shape with side length 
5 /r'Mpc for the subsequent lensing analysis. 

3. Lensing properties 

3.1. Ray-tracing 

In this section, we illustrate the techniques used to derive the 
strong lensing properties of the clusters in the MareNostrum 
Universe cosmological volume. The deflection angle maps are 
calculated as explained in several previous papers (see e.g. 
Meneghe tti et al.pOOO 2005 i. The particles in each cube are 



used to produce a three-dimensional density field, by interpo- 
lating their position on a grid of 5 12 3 cells using the Triangular 
Shaped Cloud method (Hockney & E astwood|198 8i). Then, we 
project the three-dimensional density field along the coordinate 
axes, obtaining three surface density maps used as lens 
planes in the following lensing simulations. 

The following step consists of tracing bundles of light rays 
through a regular grid covering the central part of each lens 
plane. We choose to set the size of this region as 1 .5 x 1 .5 hr 1 
Mpc 2 comoving. This choice is driven by the necessity to study 
in detail the central region of the clusters, where critical curves 
form. However, we do this by taking into account the contri- 
bution from the surrounding mass distribution to the deflection 
angle of each ray. 

We first define a grid of 256 x 256 "test" rays, for each of 
which the deflection angle is calculated by directly summing 
the contributions from all cells on the surface density map 



\Xh,k - Xi 



where A is the area of one pixel on the surface density map 
and Xh.k and jcy are the positions on the lens plane of the "test" 
ray (h, k) and of the surface density element (i, J). Following 
Wambsganss et al. ( 1998}, we avoid the divergence when the 



distance between a light ray and the density grid-point is zero 
by shifting the "test" ray grid by half-cells in both directions 
with respect to the grid on which the surface density is given. 
We then define an higher resolution grid of rays covering the 
same region. We determine the deflection angle of each new 
ray by bi-cubic interpolation between the four nearest test rays. 
The grid size is chosen such that the resolution of the deflection 
angle map is fixed at 0.2". This results into ~ 2400 x 2400 rays 
traced through the central region of deflectors at z\ - 0.3. The 
number of grid point obviously increases at lower redshift due 
to the combined increase of the physical size of the grid and 
to the decrement of the angular diameter distance. Conversely, 
the number of pixels decreases at higher redshifts. We selected 
an upper bound of 3500 x 3500 grid points because the time 
consumption and the needed memory for the code become too 
demanding, and a lower bound of 650x650 grid points in order 
to be able to fairly capture the relevant structures in the deflec- 
tor. 

3.2. Strong lensing clusters 

For a fixed source redshift, a cluster can produce strong lensing 
events if it develops critical lines on the lens plane. These lines 
correspond to the caustics on the source plane. Only sources 
within the caustics have multiple images. Only sources which 
happen to lay close to the caustics are strongly distorted and 
magnified. 

A cluster can produce both tangential and radial critical 
lines, where the tangential and the radial magnifications di- 
verge, respectively. The critical lines form where 



k(x) + y{x) = 1 



(2) 



where k(x) and y(x) = [71,72] are the convergence and the 
shear at the the position x on the lens plane. Both the con- 
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vergence and the shear are linear combinations of the spatial 
derivatives of the deflection angle components, a - [01,02]' 



(3) 
(4) 
(5) 



It can be shown that the convergence is the surface density di- 
vided by a critical surface density: 
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(6) 



which depends on the angular diameter distances between ob- 
server and lens, D\, between the lens and the source, D\ & , and 
between the observer and the source, D s : 



^-crit — 



4nG D ts D s 



(7) 



It is clear that, in order to be a strong lens, the cluster conver- 
gence and shear must be large enough, such that their sum is 
larger than unity somewhere. 

In order to focus on the subsample of strong lenses, we start 
by selecting those halos that are capable of developing critical 
lines. For doing this, we use the 256 x 256 test rays first. We 
numerically calculate the spatial derivatives of the deflection 
angles to compute k and y and to determine the positions of the 
critical points. If at least a critical point is found in the three 
cluster projections using these maps, we consider the cluster 
for further, more detailed analysis. We are aware that, using this 
selection criterium, all lenses whose critical lines have sizes 
smaller than ~ 20/T 1 kpc comoving, corresponding to the spa- 
tial resolution of the coarse grid, are not included in our analy- 
sis. Because of this limitations our results should be used with 
caution in reference to small strong lensing systems. Instead, 
such objects would never produce significant image splittings 
and large distortions or giant arcs, which are the most relevant 
strong lensing features in this work. 

As a result of this preliminar selection, it turns out that 
49, 366 clusters produce critical curves for sources at redshift 
z s = 2 in at least one of their projections. For these objects, 
we repeat the calculation of the deflection angles on grids with 
higher spatial resolution on all the three projections, obtaining 
148, 098 deflection angle maps. 

3.3. Cross-sections for giant arcs 

Once the high-resolution deflection angle maps for the three 
projections of each numerical cluster are computed, the strong 
lensing efficiency for long and thin arcs was evaluated by us- 
ing the fast, semi-analytyc algorithm presented in |Fedeli et al.| 
2006). The reader is referred to the quoted paper for details, 



while here we give just a quick overview of the method. The 
lensing efficiency is quantified by the lensing cross section for 
highly distorted arcs o~d . This is defined as the area of the re- 
gion surrounding the caustics within which a source is mapped 



on the lens plane as an image with a minimal length-to-width 
ratio d(,. The size of the lensing cross section is related to the 
expected number of arcs with a minimal distortion observed 
behind the cluster. In fact, the number of arcs above a mini- 
mal surface brightness So expected from a cluster with cross 
section cr for sources at redshift z s is 



IT 



o-{z s )n s (S,Zs)dSdzs 



(8) 



where n s (S,z s ) is the number density of sources with surface- 
brightness S and redshift z s . Note that lensing does not change 
the source surface brightness, thus N arcs does not depend on 
the magnification, at least if the PSF size is smaller that the 
size of arcs. This generally applies for extended arcs like those 
considered here. 

When sources are much smaller than the characteristic 
length over which the lensing properties of the deflector change 
significantly, they can be considered as pointlike. In this case 
the lens mapping can be linearised and the length-to-width ra- 
tio of the distorted images, d, is simply given by the ratio of the 
eigenvalues of the Jacobian matrix at image position. Hence, 
in this case the cross section for arcs with length-to-width ratio 
larger than some threshold do is by definition the area of the 
lens plane where the eigenvalue ratio is larger than do, mapped 
back to the source plane. This framework can be also easily 
modified to account for the extended size of real sources, by 
convolving the lensing properties over the typical source do- 
main, assumed here to be circular with angular radius of 0.5" . 
Several studies have shown that the properties of the sources 
are relevant for determining the shape of gravitational arcs (see 
e.g. |Meneghetti et al.|20"08"t|Gao et al.|2009| l. Our method can- 
not take into account all the effects due to source morphologies 
and luminosity profiles, however the intrinsic ellipticity of real 
sources is accounted for according to the elegant algorithm pro- 



posed by Keeton ( 2001 1. 

Because of the huge number of cross sections computed in 
this work, we focused only on a single value for the length- 
to- width threshold, namely do = 7.5. While a distribution of 
thresholds would be preferable, we expect the change in do to 
produce only a shift in the normalisation of cross sections, by 
leaving every qualitative conclusion unchanged. As mentioned 
above, we consider one source redshift, z s = 2 only. 

We computed the strong lensing cross sections for each of 
the 49, 366 high resolution deflection angle maps, produced as 
described in section |3.1| A large number of these however is 
vanishing, due to the fact that, even though the deflector is able 
to produce critical curves, these are small compared to the typ- 
ical source size to efficiently distort images. It turns out that 
only 6375 clusters have at least one projection with non van- 
ishing cross section for z s = 2. In total, the cluster projections 
capable of large distortions are 11,347. The lensing cross sec- 
tions range from a minimal value of 9.5 x 10~ 8 /z~ 2 Mpc 2 to a 
maximal value of 1.9 x 10~ 2 /T 2 Mpc 2 (the properties of this 
super-lens are shown in Appendix K). However, the vast ma- 
jority of the lenses capable of producing giant arcs have cross 
sections larger than ~ 10~ 4 /i~ 2 Mpc 2 , as it can be seen in Fig. [I] 
where we show the distribution of the lensing cross sections for 
giant arcs among the clusters analyzed here. 
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Fig. 1. The distribution of the lensing cross sections for giant 
arcs of all the strong lensing clusters between z\ — and Z\ = 2 

Given the large number of clusters analyzed here and the 
computational time required to analyze them it was not pos- 
sible to calculate the lensing cross sections for several source 
redshifts. Thus it is not easy to convert the lensing cross sec- 
tion into a number of arcs using the Eq. [8] Nevertheless, we 
can estimate such number using some approximation. If we 
assume that the lensing cross section evolves with redshift as 
cr(z s ) = cr x /o-(Zs), where cr is the lensing cross section for 
sources at z s — 2 and = cr(z s )/cr is a scaling function, then 
the number of arcs detectable behind a cluster can be expressed 
as follows: 



Narcs(So) = 0"X 



XT 



crx « eff (5 ) 



fo-(Zs)n s (S,z s )dSdz s 



(9) 



In the last equation, we have introduced the effective source 
number density. Apart from the dependency on the scaling of 
the lensing cross section with the source redshift, which will 
be discussed below, the effective source number density is set 
by the minimal surface brightness (i.e. flux per square arc- 
sec) of detectable arcs. Thus, it is determined by the charac- 
teristics of the observation, i.e. by the throughput of the in- 
strument and by the level of the background. Using the op- 
tical simulator SkyLens (Meneghetti et al.||2008} Meneghetti 
et al. 2009, in prep.), we have simulated a deep exposure 
of 8000s with the Advanced Camera for Surveys on board 
the Hubble Space Telescope in the F775W filter. This code 
uses the morphologies, the luminosities, and the redshifts of 
the galaxies in the Hubble Ultra-Deep-Field ( |Beckwith et al.| 
2006) to produce extremely realistic images of the sky includ- 
ing several observational noises. Setting the background level 
to 22.4 mag arcsec~ 2 , the lcr and the 3cr detection thresholds 
above the background r.m.s. are So ~ 25.78 mag arcsec~ 2 and 
So ~ 24.58 mag arcsec~ 2 , respectively. The galaxy number 
counts per square arcmin above the detection thresholds are 
shown in the left panel of Fig.|2]as a function of redshift. Here, 
we use the photometric redshifts of the HUDF galaxies from 



the public catalog by Coe et al. (2006). 



The scaling function f a is expected to depend on several 
properties of the lenses, like their redshifts, density profiles, 
ellipticity, and substructures. Therefore, adopting a universal 
scaling law is certainly a gross approximation. On the other 
hand, it is useful also to have a rough estimate of the effective 
number density of sources in order to link a quantity like the 
cross section, which is not directly measurable, to something 
which can be observed, like the number of arcs behind a cluster. 
To estimate the scaling function f lT for a cluster at reshift z\, we 
use a toy lens model with an NFW density profile and fixed 
projected ellipticity e = 0.2. The ellipticity is introduced in 
the lensing potential as discussed in Meneghe tti et al.| ([2003b). 
Using the same algorithm used to analyze the deflection angle 
maps of numerically simulated clusters, we measure how the 
lensing cross section grows as a function of the source redshift. 
In the right panel of Fig. [2] we show the scaling functions for 
a cluster with mass 10 15 lr x M & at several redshifts between 
z = 0.2 and z = 0.8. 

The effective number counts derived as explained above 
are reported in Tab. [T] for different cluster masses and red- 
shifts. First, we note that the dependence on the mass is weak, 
which allows to extend the validity of these calculations to a 
broad range of masses. Second, the rise of the scaling function 
for increasing lens redshift compensate for the lower number 
of galaxies at high redshift. Thus the effective source number 
counts do not drop, but tend to increase as the lens redshift 
increases. Using the Eq. [9] we can finally link the number of 
arcs expected for a given lensing cross section cr to the effec- 
tive number density of background sources, i.e. to the depth 
of the observation. For example, for a cluster with cross sec- 
tion 10~ 3 lr 2 Mpc 2 the expected number of giant arcs varies 
from 0.3 to 1.6 for n e ff in the range [40 - 200]. Conversely, in 
Fig. [3] we show the lensing cross section required for A^cs = 1 
as a function of the effective number density of background 
sources. Even for very high effective number counts (or equiv- 
alently very deep exposures), the lensing cross section needs to 
be very large in order to expect at least one arc behind a galaxy 
cluster. For example, for a cross section of cr — 10~ 3 /T 2 Mpc 2 
the effective number density of background sources needs to be 
~ 130 galaxies per square arcmin. Such number density needs 
to be doubled in order to expect to observe two arcs, and so 
forth. 

In the rest of the paper, we will use the lensing cross 
section to discriminate between lenses of different strengths. 
Observationally, the lensing cross section is not a directly mea- 
surable quantity. A possible method to estimate the lensing 
cross section is through the detailed parametric reconstruction 
of the lens potential (Meneghetti et al. in prep.). Indeed, the de- 
flection field can be readily derived from the lensing potential 
and used to measure the lensing cross section using the same 
method adopted in our simulations. Other methods to define 
the strength of the lenses may be based on quantities which 
are more directly measurable, like the angular separations of 
multiple images, which can be used to estimate the size of the 
critical lines, etc. However, in this paper, given the huge size of 
the cluster sample considered, we could not explore this other 
possibility. 
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Fig. 2. Left panel: Galaxy number counts per square arcmin in redshift bins as obtained by simulating a deep observation of 
8000s with HST/ACS in the F775W filter. The two histograms refer to two different SExtractor detection thresholds, namely 
So = 25.78 and So = 24.58 mag arcsec~ 2 , which correspond to lcr and 3cr above the background r.m.s. See text for more details. 
Right panel: scaling functions for cross sections of lenses with mass M = 10 15 h~ l M e at four different redshifts: 0.2 (solid line), 
0.4 (dotted line), 0.6 (dashed line), and 0.8 (dot-dashed line). 
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Table 1. Effective galaxy number counts per square arcmin be- 
hind clusters with different masses and redshifts. The counts 
refer to an observation with HST/ACS in the F775W filter with 
an exposure time of 8000s. In each column, the biggest and 
the smallest number correspond to detections at lcr and at 3cr 
above the background r.m.s. 
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by 



is given 



Fig. 3. The lensing cross section required for an expectation 
value of one giant arc behind a cluster as a function of the ef- 
fective number of background sources. 



4. Cluster masses 

The easiest way to characterize a cluster lens is through its 
mass. In the following, we refer to the cluster mass as the mass 
contained in spheres of radius R vn . Such virial radius encloses 
a density of A v j, (z) times the closure density of the Universe at 
the redshift of the cluster, 



Pc(z) = 8^G [°'" (1+z)3 



■] 



(10) 



where Ho is the present value of the Hubble constant, z the red- 
shift, G is the gravitational constant. The virial overdensity A v i r 
depends on redshift (see Gunn & Gott|1972" Bryan & Norman 



M vil = -7tR* k p c (z)A v iAz) 



(11) 



In alternative to the virial mass, different mass definitions are 
often adopted in literature, such as the mass corresponding to a 
constant over-density A, with A = 200, 500, or 2500 times the 
critical density. The general definition in this case is 



M A = -nR 3 A p c (z)A . 



(12) 



The mass function of objects in the MareNostrum 
Universe is in very good agreement with the theoretical ex- 
pectations ( [Sheth & Tormen|2002|[Gottlober et al.|2008) >. This 
is shown in Fig. |4] where the number of halos above a minimal 
virial mass is shown for different redshifts and compared to the 
predictions of the Sheth & Tormen mass function at redshift 
Z = 0. At this redshift more than 4,000 cluster sized objects 
with masses larger than 10 14 /i _1 M o are found. About 58,000 
objects have masses larger than 10 13 h~ l M B . At redshift z - 1 
more than 30,000 objects with masses larger than 10 13 h~ l M Q 
are detected. 
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Fig. 4. The number of halos above a given mass within the 
MareNostrum Universe simulation box at four different red- 
shifts between and 2. The dotted line shows the theoretical 
expectations from the Sheth & Tormen mass function. 



Selecting the clusters via their strong lensing efficiency im- 
plies that only the high mass tail of the distribution is prop- 
erly sampled. Indeed, since the amplitude of the gravitational 
deflection directly depends on the mass, strong lenses are the 
most massive objects at each epoch. In particular, we expect 
a minimal mass below which clusters do not develop criti- 
cal lines and are unable to produce very distorted images like 
gravitational arcs. Since clusters must be located at a conve- 
nient angular diameter distance between the observer and the 
sources, the number and the typical mass of strong lensing clus- 
ters should vary as a function of both the lens and the source 
redshifts. Assuming a fixed source redshift of z s = 2, the mass 
distribution of strong lensing clusters in a comoving volume 
of 500 3 /;~ 3 Mpc 3 at different redshifts is given in the left panel 
of Fig. [5] The color levels show the number counts of critical 
clusters in the M v ; r - z plane. Lighter (darker) colors corre- 
spond to smaller (larger) number counts. The outer dotted con- 
tour show the limits of the distribution: no critical clusters have 
been found outside the region enclosed by this line. The two 
inner contours correspond to the 50% and to the 90% of the 
peak of the distribution. Thus, they show how rapidly the criti- 
cal cluster counts decrease as a function of both mass and red- 
shift. As the figure shows, the region of the plane where clus- 
ters are able to produce critical lines extends down to masses of 
groups at the most favorable redshifts. However these are very 
rare objects. At redshifts larger than 1.2, or smaller that 0.2, 
the mass threshold grows rapidly, while the number counts of 
critical clusters drop. These lenses are too close to the sources 
or to the observer to be critical. Requiring that clusters are also 



able to produce giant arcs produces an additional selection ef- 
fect. Using the same convention as for the dotted contours, the 
solid contours refer to arcs with non-vanishing cross-sections 
for giant arcs. Note that there are no clusters at zi > 1.3 that are 
able to produce large distortions, although they could be still 
efficient for sources at much larger redshift. The most massive 
clusters in the box are still able to develop small critical lines 
up to z\ = 1.7. The bulk of clusters producing giant arcs is con- 
centrated at 0.15 < z\ < 0.8. Finally, the dashed contours show 
the distribution of the lenses with lensing cross sections above 
10~ 3 /i~ 2 Mpc 2 . These lenses are likely to be the most easily 
targeted for strong lensing studies and contribute significantly 
to the lensing signal in the universe given that, as discussed 
in the previous section, the expected number of strong lens- 
ing features produced by these objects is by far larger than for 
lenses with smaller cross sections. These, on the other hand, are 
more abundant, thus they will dominate the total lensing optical 
depth. As shown by the contours, these objects are confined in 
a smaller area on the M v ; r - z plane. They are typically clusters 
with masses exceeding few times 10 14 h~ l M B and with redshift 
below unity. Most of them are concentrated in a narrow redshift 
window between 0.2 < z < 0.6. 

Although we see an interesting selection effect in the 3D 
masses, what really matters for strong lensing is the projected 
mass. In particular, the mass contained in a cylinder around 
the cluster center, where the critical lines form. For each clus- 
ter in our sample, we measure the projected mass within /?2500, 
for each of the three projections used for ray-tracing. 
This is defined as that of a sphere encompassing a mean den- 
sity of 2500 x p c . Typically, it corresponds to a region which 
is large enough to contain the cluster critical lines. Here, the 
projected mass is obtained by integrating all the mass in a 
cylinder of height 5/r'Mpc. We show in the right panel of 
Fig.[5]the distribution of the strong lensing clusters in the plane 
_ Z- Interestingly, although it appears clear that strong 
lensing depends on the mass in the central region of the de- 
flectors, the spread in projected mass is larger by about one 
order-of-magnitude than in 3D. There are clusters which have 
relatively small mass projected in the core but that are still ca- 
pable to produce strong lensing effects of different intensity. We 
interpret this result as due to the importance that other proper- 
ties of the lenses have for strong lensing, like the amount of 
substructures and the level of asymmetry and ellipticity in the 
cluster cores, as shown in Meneghetti et al. (2007a). In several 
cases, and especially for clusters producing mild strong lens- 
ing effects (i.e. clusters with critical curves) the excess of shear 
produced by a clumpy and asymmetric mass distribution can 
compensate for the low value of the central convergence. 

5. Halo triaxiality and orientation 



Since simulated galaxy clusters are triaxial (see e.g. Jing & 



Suto|2000| |Gottlober & Yepes|2007| , their projected mass de- 
pends on their orientation with respect to the line of sight. In 
order to evaluate how this impacts on the strong lensing ability 
of clusters, we measure their triaxial best fit model and discuss 
the correlation of the orientation with the occurrence of critical 
lines and giant arcs. To do that, we measure the moment of iner- 
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Fig. 5. The distributions of strong lensing clusters in the M v ; r - z (left panel) and in the M p 25oo - z plane (right panel). The 
color levels correspond to different number counts of halos in the MareNostrum Universe (500 3 /;~ 3 Mpc 3 comoving) which are 
critical for a source redshift of z s = 2. The inner dotted contours show the levels corresponding to 90% and 50% of the peak of 
the distribution. The outer dotted contour encloses 100% of the halos, thus it shows the minimal and the maximal mass of halos 
producing critical lines for each redshift bin. Similarly, the solid and the dashed contours refer to the distributions of the halos 
with cross sections for giant arcs larger than and than 10~ 3 /T 2 Mpc 2 . The color-bar on the top of each panel show the link 
between the colors and the log of the number of cluster per mass and redshift bin. 



tia tensor In of each cluster in the sample. The cluster particles 
within R v i r are sorted in a regular cubic grid of 256 x 256 x 256 
cells. The mass density in each grid cell is then computed The 
inertial tensor components are given by 

Iij = Y j m k (RlS i j-R Ki R kJ ), (13) 

k 

where m k is the mass of in the fc-th selected cell and R = [Ri] is 
the vector which identifies the cell position with respect to the 
centre of mass of the system. The triaxial model of the cluster 
and its principal axes are obtained by diagonalizing the inertial 
tensor, finding its eigenvalues and eigenvectors. 

The fits show that clusters have prolate triaxial halos and 
the distributions of the axis ratios of strong lensing clusters is 
not significantly different from that expected for the general 
cluster population (see e.g. Jing & Suto 2002). This is in agree- 
ment with the results of Henna wi et al.| (|2007 ), who also find 
that strong lensing clusters are not significantly more triaxial 
than normal clusters. However, strong lensing clusters seem 
to be affected by an orientation bias. In Fig. [6| we show the 
cumulative probability distribution function of the angle be- 
tween the major axes of the inertial ellipsoid and the line of 
sight to the cluster. The solid and the dot-dashed lines indi- 
cate the results for the whole sample of lensing clusters and for 
the sub-sample of clusters producing giant arcs. The dashed 
line shows the distribution of the orientation angles of the most 
efficient lenses in the sample, i.e. with lensing cross section 
cr > 10~ 3 /;~ 2 Mpc 2 . We also display, using a dotted line, the 
distribution corresponding to totally randomly oriented lenses. 
We find that lensing clusters tend to be preferentially aligned 
with the line of sight. This orientation bias increases with the 
strength of the lens: the median angle for critical clusters is 
~ 57 deg, while for the sub-sample of clusters capable of pro- 
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Fig. 6. The cumulative probability distribution function of the 
angles between the major axes of the strong lenses in the 
MareNostrum Universe and the line of sight. Shown are the 
results for the clusters with critical lines, (solid line), for the 
clusters with cross section for giant arcs larger than zero (dot- 
dashed line), and for the clusters with cross section for giant 
arcs larger than 10~ 3 /T 2 Mpc 2 (dashed line). The dotted line 
shows the expected distribution in case of randomly oriented 
halos. 

ducing giant arcs the median angle is ~ 50 deg. The median 
decreases to ~ 47 deg for clusters with lensing cross section 
cr > 10~ 3 /T 2 Mpc 2 . Note that in the case of random orientation 
we should expect a median angle of 60 deg. This is an important 
effect, which can affect the conclusions of many studies aim- 
ing at estimating the mass of clusters through strong lensing or 
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at measuring cosmological parameters using the abundance of 
highly elongated arcs on the sky. In fact, we expect that, due 
to the orientation bias, 3D strong-lensing masses are biased 
high, if the approximation of spherical symmetry is used to 
convert the measured two-dimensional into three-dimensional 
mass profiles. Moreover, this alignment bias has to be properly 
modeled when estimating the lensing optical depth for a pop- 
ulation of strong lenses in a given cosmology. Similar results 



were found by Hennawi et al. (2007 1. They also find a corre 



lation between strong lensing and orientation of the lenses and 
they find similar distributions of the orientation angles as those 
we find here. 

Apart from the orientation, the halo triaxiality is important 
because it determines the projected shape of the lenses. It has 
been shown in several papers that, for a fixed mass, the strong 
lensing cross section is larger for higher ellipticities of the pro- 
jected mass distribution (Menegh etti et al.|2003b| |MeneghettT| 
et al.|2007a] l. With such a large sample of lensing clusters, we 



can address the question of what is the distribution of their pro- 
jected ellipticities. These are measured similarly to the three- 
dimensional shape of the lenses. We measure and diagonalize 
the inertial tensor of the cluster mass distribution projected on 
a regular grid of 256 x 256 cells. We select those cells where 
the surface density exceeds some thresholds. The thresholds 
we use are given by the mean surface densities at R Y [ r and at 
0.1 xR v i r . Thus, we measure the projected ellipticity both in the 
outer and in the inner cluster regions. The probability distribu- 
tion functions of the projected ellipticity are shown in Fig. [7] 
The ellipticity is defined as e = (1 - b/a)/( 1 + b/a), where a 
and b are the major and minor axes of the ellipse. The solid, 
the dot-dashed, and the dashed lines refer to critical clusters 
(assuming again a source redshift of z s = 2), to clusters with 
non-vanishing cross section for giant arcs, and to clusters with 
large cross section for giant arcs (> 10~ 3 /T 2 Mpc 2 ), respec- 
tively. The left and the right panels show the distributions of 
the outer and of the inner ellipticities. We find that the pro- 
jected cores are more elliptical, with distributions which peak 
at e ~ 0.2 - 0.4. It is interesting to note that critical clusters 
have a bimodal ellipticity distribution: a large number of clus- 
ters have extremely elongated cores with ellipticities which ex- 
tend to e = 0.9. Since we are fitting each lens with a single 
ellipse, these are mainly clusters with substructures near the 
centers which mimic large ellipticities. We shall recall that the 
tangential critical lines form where 

K + y=l (14) 

As discussed in Torri et al. ( 2004} , the shear produced by the 
substructures enhance the the ability of the clusters to pro- 
duce strong lensing, because it makes critical even those lenses 
where the convergence is not large enough to ensure it (k > 1 
at some point). However, although the additional shear allows 
several clusters to have critical lines, several of them still are 
unable to produce large distortions. Indeed, the second peak 
in the PDFs at large ellipticities is less prominent for clusters 
producing giant arcs, and even less for clusters with large lens- 
ing cross sections. For this class of lenses, it is more impor- 
tant to have a large amount of mass projected onto cluster core. 
Thus, either they have substructures closer to center, or they 



are more strongly aligned with the line of sight, thus appearing 
a bit rounder than other critical lenses. Note that, as long as 
the lensing cross section grows, the ellipticity distribution be- 
comes unimodal but its peak shifts towards larger ellipticities, 
also suggesting that most of these cluster have substructures 
close to their centers. 

At large radii, we find that the ellipticities are smaller and 
the projected ellipticity becomes smaller as the strength of the 
lens increases. This is clearly related to the orientation bias 
discussed above. The strongest lenses typically are elongated 
along the line of sight, thus they appear rounder on the sky. 

6. Concentrations 

Several previous studies have discussed the importance of the 



halo concentration for lensing. Using simulations, Hennawi 
|et al.| ( [2007] l find that concentrations of lensing clusters are on 
average ~ 18% larger than the typical clusters in the universe. 



Broadhurst et al. ( 2008 ) report very high level of mass concen- 



trations (c ~ 10) in a sample of four well known strong lensing 



clusters. Fedeli et al. ( 2007b I show that, at a given mass, the 
strong lenses are ~ 10% to ~ 25% more concentrated than the 
average. Here, we discuss the concentrations of the clusters in 
the MareNostrum Universe. 

We measure the concentrations by fitting the density pro- 
files of the clusters in our sample with the Navarro-Frenk- 
White ( |Navarro et al.|1997| > formula, 

Ps 



P(r) = 



r/r s (l - r/r s ) 2 



(15) 



where p s is a characteristic density and r s is the scale radius. 
The concentration is defined as c — R2oo/r s . 

Instead of fitting the density profiles of individual halos, 
which are noisy, we prefer to fit the stacked profiles of clusters 
with similar redshifts and masses. The mass bins are equally 
spaced in logarithmic scale. We stack the profiles of all clusters 
in the mass bins and perform the NFW fit. Again, we select 
those objects which exhibit critical lines for sources at z s = 2, 
and, among them, those halos which are also able to produce 
large distortions, with lensing cross section above some mini- 
mal value. The concentrations of the clusters in these two sub- 
samples are compared with those of general clusters, regard- 
less of their ability to behave as strong lenses. The results are 
shown in Fig. [8] The left and the right panel refer to critical 
and to large distortion clusters, respectively. The color inten- 
sity reflects the amplitude of the concentrations. The concen- 
trations are normalized to those of general clusters of similar 
mass and redshifts. The labels in the overlaid contours indicate 
the numerical value of the normalized concentration at the cor- 
responding color level. Strong lensing clusters at moderate red- 
shifts have concentrations similar to those of the general clus- 
ter population. Only small mass lenses have a relatively small 
concentration bias (< 20%). The bias become more signifi- 
cant at small and high redshifts, where it affects also the largest 
masses. Due to their short distance to the observer or to the 
sources, these clusters need to be very concentrated in order to 
focus the light from distant sources. Note that the bias is mass 
dependent. As the mass decreases, the bias is stronger. This is 



10 



M. Meneghetti et al.: Strong lensing by clusters in the MareNostrum Universe 




Fig. 7. The probability density functions of the projected ellipticity of the lensing clusters in the MareNostrum Universe. Shown 
are the results for the clusters with critical lines (solid lines), for the sub-sample of clusters with non-vanishing cross sections for 
giant arcs (dot-dashed lines), and for cluster with cr > 10" 3 /T 2 Mpc 2 . The left and the right panels refer to ellipticities measured 
at /?vi r and at 0.1 x R vir . 




Fig. 8. NFW concentrations as a function of the halo mass and redshift. Left and right panels refer to clusters with critical lines 
and to clusters with non-vanishing lensing cross section for giant arcs, respectively. The concentrations are normalized to those 
measured on the whole sample of clusters in the simulation box, regardless of their ability to produce strong lensing effects. 
Different colors are used to encode the different values of the normalized concentrations. The labelled contours provide the link 
between the color scale and the concentration values. 



a clear selection effect: if we require a cluster to be critical or 
even to produce large arcs, only the most concentrated halos in 
the smallest mass bins are able to satisfy the requirement. 

As mentioned above, lensing probes the projected mass 
distribution of clusters. The concentrations are typically mea- 
sured by fitting multiple image systems and arcs with com- 
binations of projected parametric models. Then, the three di- 
mensional density profiles are determined by assuming spher- 



ical symmetry. As we discussed earlier, clusters have triaxial 
shapes, thus the assumption of spherical symmetry is generally 
wrong. Moreover, as we have shown in the previous section, 
strong lensing clusters tend to be seen along their major axes. 
For these reasons, the concentrations measured in 2D through 
strong lensing are expected to be more biased compared to 3D 
concentrations. This effect is also discussed in Hen nawi et al.l 



(2007), where a comparison of 2D vs 3D concentrations of in- 
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Fig. 9. Same as in Fig. [8| but showing 2D-concentrations of lensing clusters. The 2D-concentrations have been normalized to the 
3D-concentrations of the whole cluster sample. Starting from the top-left panel, we show the 2D-concentrations as a function of 
mass and redshift for 1) clusters with critical lines; 2) clusters with non-vanishing lensing cross sections for giant arcs; 3) clusters 
with cross section cr > 1CT 3 /i~ 2 Mpc 2 ; 4) clusters with cross section <x > 2 x 1CT 3 /j~ 2 Mpc 2 . 



dividual clusters led to the conclusion that the former are typ- 
ically ~ 20% larger than the latter. To verify this result, we 
proceed to fit the surface density profiles of our strong lensing 
clusters in their projections. Again, to do this, we stack the pro- 
files in mass and redshift bins. The fitting formula is given by 



the truncated NFW surface density profile (Meneghetti et al. 
[2000) 1, 



r 

Jo 



,(r)df , 



(16) 



where £ is the coordinate along the line of sight and £ is the 
component of r perpendicular to The maximum of £ is given 
by half the size of the sub-box containing each cluster, i.e. 
2.5 /r'Mpc comoving. Using the dimensionless coordinate 



on the projection plane x = £/r, and defining the quantities 
u = arcsinh(£/£), and k s = p^E" 1 , the previous equation can 
be written as 



^nfw( x ) 
where 

fix) = - 



2K s £ cr /0) 



(x 2 - 
1 



arctan 



1)3/2 

x sinh u 



x- 1 



Vx 2 - 1 



: tanh 



x 2 - 1 1 + x cosh u 



if x > 1 ; 



/(*) = 



2cosh(|)sinh(f) 4cosh(|) 3 sinh(f) 
3(l+cosh«) 2 3(l+cosh«) 2 



(17) 



(18) 



(19) 
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if x — 1 ; and 



-arctanh 



(l-*2) 3/2 

1 x sinh u 



1 



Vl -x 2 



: tanh 



(I) 



1 - x 2 1 + x cosh u 



(20) 



if x < 1. In the previous formulae, M max = arcsinh(^ max /^). 

The resulting 2D concentrations, for several classes of 
strong lenses, as a function of mass and redshift are shown 
in Fig. [9] As done in Fig. [8] the 2D-concentrations are nor- 
malized to the 3D-concentrations of general clusters of simi- 
lar masses and redshifts. critical and large distortion clusters 
separately and we normalize the measured concentrations with 
those derived fitting the 3D density profiles of the whole sam- 
ple of clusters in the cosmological box. Starting from the top- 
left panel we show the results for clusters with critical lines 
and for clusters with lensing cross sections for giant arcs larger 
than 0, 1(T 3 /T 2 Mpc 2 , and 2 x 1CT 3 /T 2 Mpc 2 , respectively. 
As expected, the bias grows compared to the 3D case, and the 
amount by which it increases depends on the class of lensing 
clusters we are considering. For critical clusters, the ratios be- 
tween 2D concentrations and the corresponding 3D concentra- 
tions are of the order of ~ 1.2 for intermediate redshift clus- 
ters but they become larger than 1.3 at low and at high red- 
shifts. For clusters which are able to produce giant arcs, the 2D- 
concentration bias is significantly larger. As discussed for the 
3D concentrations the amplitude of the bias depends on both 
redshift and mass: smaller masses at short distances from the 
observer or from the sources have the largest biases. Moreover, 
increasing the lensing cross section for giant arcs, the concen- 
tration bias grows dramatically. For example, massive clusters 
(M ~ 10 15 /z _1 M ) at the most efficient redshifts for strongly 
lensing sources at z s = 2 (z\ ~ 0.4) with lensing cross sec- 
tions larger than 2 x 10~ 3 /i~ 2 Mpc 2 have 2D-concentrations 
which are typically larger by ~ 50% than the 3D-concentration 
of general clusters. For smaller masses and redshifts, the 2D- 
concentrations can be higher than expected in 3D by more than 
a factor of two. This is a consequence of the orientation bias 
discussed in the previous section. In order to be able to pro- 
duce large and very elongated arcs, clusters laying too close to 
the observer or to the source must be optimally oriented and 
extremely concentrated. Due to triaxiality, the concentrations 
measured from the 2D mass distributions of these clusters are 
much bigger than the corresponent 3D concentrations (Oguri 
|et al.|2005||Gavazzi|2005l l. Note that, as discussed in Sect. [33] 
a lensing cross section of cr ~ 10~ 3 /i~ 2 Mpc 2 corresponds to 
an expectation value of ~ 1 giant arc in a deep HST observa- 
tion. Clusters like A1689, which contains about 10 arcs with 
length-to-width ratio larger than 7.5 ( Sand et al. 2005} , are thus 
expected to have extremely large lensing cross sections. If the 
properties of real clusters are well reproduced by the clusters in 
our simulations, these very efficient strong lenses are likely to 
have extremely biased 2D-concentrations, as recently reported 
by |Broadhurst et al.| ( |2008] l (see also |Oguri etliij|2i)()9"] >. Our 
findings are in agreement with the results recently published 
by Oguri & Blandford (2009 ), who use semi-analytic models of 
triaxial halos to estimate that the projected mass distributions 



of strong lensing clusters have ~ 40 - 60% larger concentra- 
tions compared with typical clusters with similar redshifts and 
masses (see also |Sereno et al.|2 010). 

7. X-ray luminosities 

Gas physics is known to be potentially very important for 
strong lensing (see e.g. Puchwein et al.|[2005| |Hilbert et al. 
2008). Several processes taking place in the Intra-Cluster- 
Medium (ICM), like cooling, heating, energy feedback from 
AGNs and supemovae, thermal conduction, etc. can also af- 
fect the distribution of the dark-matter in clusters, influencing 



the shape of the density profiles (Dolag et al. 2004; Puchwein 



|& Hilber t 2009 Yepes et al.||2007[ ) as well as the triaxiality 
of the dark matter halos. Phenomena like cooling, star forma- 
tion, and energy feedback change the thermal properties of the 
ICM, thus influencing the X-ray emissivity. In fact, several nu- 
merical studies report that X-ray luminosities in non-radiative 
simulations are higher than in simulations where cooling and 
feedback are active, while the Lx - M relation derived from 
the same simulations is less steep than observed (e.g 



Short 



|et al.|2010||Mantz et al.|2008) . The extreme complexity of the 
processes involved presents a serious challenge for simulating 



them accurately in a cosmological setting (e.g. |Borgani et al. 
2004). Nevertheless, also a non-radiative simulation like the 
MareNostrum Universe can provide useful qualitative infor- 
mation on the possible correlation between strong lensing and 
X-ray emission by galaxy clusters. Here, we focus in particular 
on the X-ray luminosity, which is often used to select clusters 
for strong lensing surveys (Luppi no et al.|19 99). 

The X-ray bolometric luminosity is calculated from the 
temperature and internal energy of each gas particle in the 
simulated clusters. In short, the X-ray luminosity is the sum 
of the contributions to the emissivity from each gas particle, 
Lx = s i> where the sum extends over all the particles within 
R v i r . The emissivity of each gas element can be written as 



e,- = n ei {n H ,iMTi,Z,)dVi 



(21) 



where and are the number densities of electrons and 
of hydrogen atoms, respectively, associated to the i-th gas par- 
ticle of given density p,-, temperature T,, and metallicity Z,. 
The cooling function A(T,Z) is calculated by using a MEKAL 
plasma model ( |Mewe et al.|[T985l [1986] |Liedahl et al"7|[T995> 
implemented in the XSPEC software package (Arnaud[1996). 
For the metallicity, we adopt the typical value Z = 0.3Z o 
( |Fukazawa e"TaT1[T998l |ScMndler|[T99"9l see also |Bartelmann| 
& White |2003] >. Finally, dVi = nii/pi is the volume of the i-th 



gas particle of mass m,. 

The distribution of strong lensing clusters in the Lx~z plane 
is shown in Fig. 10 where we use the same notation as in Fig. [5] 
Again the counts correspond to the comoving volume of the 
MareNostrum Universe. Not surprisingly, given that the X-ray 
luminosity scales with the cluster mass (e.g Kaiser 1986), the 
distribution of the strong lenses in the Lx — z plane is very sim- 
ilar to that in the M — z plane. As found for the masses, at each 
redshift there exists a minimal X-ray luminosity below which 
no critical lenses are found. The "critical" X-ray luminosity 
reaches a minimum between z ~ 0.3 and z ~ 0.5. Increasing 
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Fig. 11. The relation between the X-ray bolometric luminosity and the cluster mass. Results are shown in four different redshift 
bins, as indicated on the top of each panel. Black, red, and blue data points (and errorbars) refer to the sub-samples of critical 
and large distortion lenses. 



the minimal lensing cross section the distributions of clusters 
producing giant arcs moves upwards and shrink along the red- 
shift axis. 

We now explore more in details the scaling of the X-ray 



luminosity with the mass for strong lensing clusters. In Fig. 1 1 
we show the Lx - M relation for four redshift bins, namely 
< z < 0.25, 0.25 < z < 0.5, 0.5 < z < 0.75, and 0.75 < z < 
1. Numerical simulations are known to be poor at describing 
the X-ray properties of the cosmic structures on the scales of 



groups (Borgani et al. 2008} . Thus, we limit our analysis to 
clusters of masses M vir > l0 u h~ l M Q . The L x - M relations 
found for clusters with critical lines, with non-vanishing cross 
sections for giant arcs, and with lensing cross sections <x > 
10~ 3 /i~ 2 Mpc 2 are shown in black, red, and blue, respectively. 
We find that, increasing the strong lensing efficiency, the slope 
of the Lx- M relation changes, becoming smaller especially at 
the lowest masses. This effect is also redshift dependent, being 
more extreme in the lowest and in the highest redshift bins. It 
shows that at the least favorable redshifts for strong lensing, 
the X-ray luminosities of the strong lensing clusters tend to be 
larger than for the general cluster population, especially if the 



lens mass is relatively small. This seems to suggest that some 
cluster property rather than the mass plays an important role 
for boosting the lensing cross sections of these small lenses, 
which also influences their X-ray luminosity. As discussed in 



Torri et al. ( 2004 1 and in Fedeli & Bartelmann ( 2007 1 mergers 
are likely to be a explanation for the effect we observe here. 



8. Cluster dynamical state 

In this section, we use indicators of the virial and of the hy- 
drostatic equilibria in clusters to investigate wether there is a 
correlation between strong lensing and dynamical activity in 
the lenses. 



8.1. Virial equilibrium 

The virial equilibrium is perhaps the most natural choice when 
trying to quantify the dynamical state of a bound structure. In 
the case of self-gravitating systems, this is quantified by the 
competition between the total potential energy and (twice) the 
internal kinetic energy. 
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Fig. 10. The distribution of clusters with critical lines (dotted 
contours and color levels), with non-vanishing cross section for 
giant arcs, and with large lensing cross sections for giant arcs 
(solid contours), <x > 10~ 3 h~ 2 Mpc 2 (dashed contours) in the 
Lx - z plane. The contour description is the same as in Fig. [5] 

In agreement with the scalar virial theorem, and following 
Shaw et al. ( 2006| l, we introduce the parameter f3 for the simu- 



lated clusters, defined as 



/?=! + 



2T-S 



U 



1 



2T-S 
\U\ ' 



(22) 



where T is the internal kinetic energy, U the potential energy, 
and S a surface pressure term that arises from considering the 
structure as contained in a limited region (in this case the finite 
sphere). If A^ v ; r is the total number of particles contained in the 
sphere of radius R v [ r , then the kinetic energy is evaluated as 



1 i=i 



and the potential energy is given by 



G v tnitrij 



2 



(23) 



(24) 



In Eqs. (23 i and (24i m,, r, and v,- are the mass, the position 
with respect to the center of mass and velocity with respect to 
the center of motion of the i— th particle, respectively. 

The surface term, S , is computed by selecting the outermost 
20% of the particles that are inside R v [ r . Let us call (2o.8 this set 
of particles. If /?o.8 and Ro g are respectively the innermost and 
the median distances from the cluster center of the particles 
within Qo.s, then 



S = 



R 0.9 



vir 0.8 lEgojg 



2 m,||v,|| 



(25) 



where the sum is now extended to all the particles between 7?o.8 
and R wn . 

For an ideal cluster in perfect virial equilibrium, (3 — 0. 
On the other hand, whenever the cluster is dynamically active, 




10' 

M v!r [h"' M s ] 



Fig. 12. Top panel: the median value and the 25% and 75% 
percentiles of the distributions of the parameter (3 as a function 
of redshift. Solid, dotted, and dashed line indicate the results for 
clusters with critical lines, with non-vanishing cross section for 
giant arcs, and with cross section cr > 10~ 3 /z~ 2 Mpc 2 . Bottom 
panel: the median value and the 25% and 75% percentiles of 
the distributions of the parameter [5 as a function of the cluster 
mass. 



like accreting mass or merging with substructures, j3 < 0, be- 
cause the dynamics of the structure is dominated by the ki- 
netic energy. According to the hierachical structure formation, 
we expect that structures at high redshift have /? on average 
more negative than at low redshift, the lowest mass objects hav- 
ing the smallest /3. At the typical redshift for lensing clusters 
(z ~ 0.3 - 0.5) the mass range considered here is well beyond 
the characteristic collapsing mass (AT. = few x lO l2 /r'M ). 
Thus, all our clusters are in the process of formation and have 
negative /3 parameters. As expected, we find a significant evo- 



lution of p with redshift. In the top panel of Fig. 12 we show 



the median (together with 25% and 75% percentiles) values 
of /J as a function of redshift. Different line styles indicate the 
results for clusters with critical lines, for clusters which have 
non-vanishing cross section for giant arcs, and for clusters with 
large lensing cross sections (cr > 10~ 3 /T 2 Mpc 2 , as usual). The 
/J parameter is close to zero at low redshift, and it becomes 
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increasingly more negative going to high redshift. In particu- 
lar, the value of /3 decreases from ~ to ~ -1 from z = 0.1 
to z - 1. Note that the most efficient lensing clusters, i.e. the 
clusters with large lensing cross sections, are characterized by 
smaller ji parameters. The bottom panel shows the median and 
the 25% and 75% percentiles of jS as a function of the cluster 
mass. The differences between the cluster sub-samples become 
more pronounced at the small masses. This indicates that the 
least massive clusters with large lensing cross sections tend to 
be farther from the virial equilibrium than clusters with small 
lensing cross sections. These dynamically active clusters also 
have higher X-ray luminosity, explaining the flattening of the 



L x - M relation shown in Fig. 1 1 towards the small masses. 



8.2. Hydrostatic equilibrium 

A different method for quantifying the equilibrium state of a 
cluster consists of determining the hydrostatic equilibrium of 
the gas filling its potential well. If the system is in equilibrium 
the pressure and the gravitational forces at a given position are 
exactly counter-acting. 

Given the hydrostatic equilibrium equation, V<I> = 
— VPg/pg, we define the hydrostatic equilibrium parameter Y 
at a given distance from the cluster center as 



rw = l - w(r) 

AnGM(rY 

where M(r) is the total mass in a sphere of radius § r , and 



W(r) 



-f 

Jas, 



VP 



Ps 



dA. 



(26) 



(27) 



Note that, while O is the total potential of the cluster, P„ and p g 
are the pressure and density of the gas component. 

If hydrostatic equilibrium holds in the gas shell of radius r, 
then F(r) = 0. Depending on whether the gas is compressed or 
is not thermalized, for example due to merging with substruc- 
tures or infall of material from the external regions, T becomes 
negative or positive. As shown by |Rasia et aT] ( |2004| l, the gas 
is generally not at rest inside the cluster potential wells. Non- 
negligible subsonic bulk motions contribute to the total pres- 
sure support of the gas. This typically leads to underestimate 
the total mass of the system under the assumption of hydro- 
static equilibrium (see e.g. Ascasibar et al.||2~003| |Rasia et al 



2006] 12008} |Nagai et al.|2007||Amegrio et al.|2009| |Meneghetti 
et al.|[2009p . The effect is the largest at the largest distances 
from the cluster center, where the infall of matter is more pro- 
nounced. We evaluate the hydrostatic equilibrium parameter 
r at two different cluster-centric distances, namely R500, and 
^2500- In Figure [T3] we show the redshift evolution of T. The 
median (with 25% and 75% percentiles of the distribution) in 
each redshift bin are shown for the subsamples of critical clus- 
ters, clusters capable of producing giant arcs, and clusters with 
lensing cross section larger than 10~ 3 /z~ 2 Mpc 2 (solid, dotted, 
and dashed lines, respectively). The curves are rather flat, dif- 
ferently from what found for the /3 parameter. The results at 
different radii are consistent with the previous finding of Rasia 



0.6 



3 0.4 



0.2 




Fig. 15. The fraction / su b of structures that, at each redshift step, 
change the value of the mass inside at least a rectangular vol- 
ume (see text for details) of at least 30% (black filled triangles) 
and 50% (red empty squares) with respect to the previous step, 
as labelled in the plot. 



center. We do not see a clear dependence of T on the lensing 
cross section. Note however that, while /3 is a global indicator 
of departure from virial equilibrium, i.e. it is sensitive to any 
excess of kinetic energy within the virial radius, F is only a lo- 
cal indicator of hydrostatic equilibrium, i.e. it can only be used 
to measure a local departure from hydrostatic equilibrium. If, 
as we believe, mergers play an important role for strong lens- 
ing, this parameter is less efficient to capture them. 

It is interesting to see where the clusters in our sample are 
located in the /3 - F plane. This is shown in Fig. 14 for T mea- 



et al. ( 2004| l: the excess of pressure support due to gas bulk mo- 
tions amounts to ~ 15-20% at/?5oo, and decreases towards the 



sured at Rsa) and for four different redshift bins. Similar re- 
sults would be seen using F25oo instead of Tsoo- The contours 
in each panel show the distributions of clusters in different sub- 
samples: critical lenses (black solid lines), clusters with non- 
vanishing cross sections for giant arcs (dotted lines), and clus- 
ters with large lensing cross sections (dashed lines). The in- 
ner and the outer contours indicate the 90% and the 50% of 
the distribution peaks. As the lens strength increases, the distri- 
butions shift along the F axes towards smaller /?s. The separa- 
tions between the distribution peaks become increasingly larger 
with redshift, indicating that the dynamical activity of clusters, 
as highlighted by the departure of these structures from virial 
equilibrium, has a strong impact on their ability to produce 
large distortions especially at high redshift. This is expected 
since the dynamical activity in clusters grows by going back in 
time. 
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Fig. 13. The median value and the 25% and 75% percentiles of the distributions of the parameter T as a function of redshift. 
Solid, dotted, and dashed line indicate the results for clusters with critical lines, with non-vanishing cross section for giant arcs, 
and with cross section <x > 1(T 3 /i~ 2 Mpc 2 . From left to right, the results are shown for F measured at R^m and 7?2500, respectively. 



8.3. Mass fluctuations in the cluster core 

In order to verify once more our interpretation of the results 
discussed above, we perform an additional test to verify that 
the clusters with increasingly negative values of ft are indeed 
clusters with infalling substructures. Considering the poten- 
tial impact of cluster mergers on strong lensing, we shall take 
into account not only those sub-structures that merge with the 
main cluster clump and start orbiting around it at small dis- 
tance from its center, but also those sub-clumps that orbit at 
large distances but which transit close to cluster core perpen- 
dicularly to the line of sight. Thus, for each cluster projection 
we reconstruct the mass accretion history by measuring the pro- 
jected mass within /?2500 at different epochs. Sudden jumps in 
the curves M Pi 25oo(z) signal the passage of some substructure 
across the cluster central region. We analyze these curves as 
follows. First, we smooth the curves by removing the peaks 
with a boxcar smoothing algorithm. This provides us a smooth 
accretion mass history M smoot h(z). Then, we compare the un- 
smoothed curve M p 250ofe) to the smoothed one, and we iden- 
tify as mergers those events where M Pi 250o(z) > Kx M smoot h(z). 
For this analysis we choose K = 1.3 and K = 1.5. 

In Fig.[l5]we show the redshift evolution of the fraction / su b 
of clusters exhibiting these temporal mass variations compared 
to the total sample. As expected, the fraction of clusters with 
large projected mass variations in the core tends to increase 
with increasing redshift. The fraction of clusters where the ex- 
cess of mass within Z?2500 is at l east 30% varies from ~ 0.1 at 
z - 0.1 to up to ~ 0.7 at z > 1. 

In Fig. 14 the distribution of the "merging" clusters in the 



P~T plane is given by the blue solid contours. Clearly, merging 
clusters are typically objects with extremely negative /3 param- 
eters. Their F parameters do not significantly differ from those 
of clusters with critical lines, confirming that this parameter 
is not a good indicator of substructure accretion. Interestingly, 
the distributions of "merging" clusters in the four redshift bins 
nicely overlay the distributions of clusters with large lens- 



ing cross sections for giant arcs. This confirms our interpre- 
tation that the most efficient strong lenses in the MareNostrum 
Universe are dynamically active clusters, with substructures lo- 
cated near the cluster core. Thanks to their contribution of the 
shear and to the overall convergence, these substructures boost 
the strong lensing cross sections. 

In Fig. [16] we compare the distributions of the X-ray lu- 
minosities of all clusters in the MareNostrum Universe (solid 
histograms) and of the sub-sample of merging clusters (dashed 
histograms). Clearly, the clusters ongoing a merging phase 
have higher X-ray luminosities compared to the general clus- 
ter population (see also |Rowley et al.]|2004 1. The differences 
between the distributions increase with redshift, indicating that 
the impact of mergers is larger at higher redshifts, where clus- 



ters have smaller masses. In Fig. 17 we show the X-ray lumi- 
nosity distributions of general and merging clusters in the red- 
shift range 0.25 < z < 0.5. In the left panel, we select the clus- 
ters with masses 10 U h~ l M Q < M vh - < 5 x VQ u hT x M B . Among 
these relatively small mass clusters, those ongoing a merging 
phase are significantly more X-ray luminous. In the right panel, 
we re-plot the distributions after selecting only the most mas- 
sive clusters (5 x l0 u hr 1 M e < M vir < 3 x lO 15 /z _1 M ). In this 
case, the X-ray luminosity distributions are almost identical, 
indicating that most of the differences between merging and 
general clusters appear at low masses. The results shown here 
help explaining the behavior of the luminosity-mass relations 



displayed in Fig. 1 1 As the lens strength increases, the strong 
lensing cluster population tends to be dominated by merging 
clusters, which are characterized by higher X-ray luminosities. 
For this reason, the L x -M relation tilts at the smallest masses, 
because clusters with high lensing efficiency and small mass 
are typically merging objects. 

9. Summary and conclusions 

In this paper, we have investigated the properties of ~ 
50000 strong lensing clusters at different redshifts in the 
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Fig. 14. The distributions of lensing clusters with increasing ability to produce large distortions in the yS-Tsoo plane. The contours 
indicate the levels corresponding to 90% and 50% of the distribution peaks. The black solid contours show the distribution of 
all clusters which have critical lines for sources at z s = 2. The blue-dotted and the red-dashed contours show the distributions 
of clusters with minimal lensing cross sections for giant arcs <x = and cr = 10~ 3 /i~ 2 Mpc 2 , respectively. Finally, the blue-solid 
contours indicate the distribution of clusters identified as "mergers", i.e. for which a variation of mass > 30% between consecutive 
snapshots has been measured in the projected region of radius /?2500 around the cluster center. 



MareNostrum Universe cosmological simulation. Projecting 
each of these clusters along three orthogonal lines of sight, we 
have considered almost 150000 lens realizations in total. With 
such a big number of objects, we can statistically characterize 
the population of strong lensing clusters much better than it 
was done in the past. Moreover, the MareNostrum Universe 
includes gas, thus it allows to investigate the correlation be- 
tween strong lensing and X-ray observables. We have classified 
the strong lenses into two categories, namely clusters which 
are critical for sources at z s = 2 and clusters which can induce 
large distortions in the images of these sources, i.e. form gi- 
ant arcs. We have explored several structural properties of the 
strong lensing clusters, like the masse, the shapes, the concen- 
trations, the X-ray luminosity, and the dynamical activity. 
Our main results can be summarized as follows: 

- strong lensing clusters are typically massive objects. Then- 
masses vary over two orders of magnitude. The minimal 
mass for strong lensing depends on both the redshift of the 



lenses and of the sources. For sources at z s = 2, we find 
that clusters can develop critical lines down to masses of 
~ 10 13 ft M Q . Requiring that clusters are also able to 
produce giant arcs, increases the mass limit by almost one 
order of magnitude. The lensing cross section further de- 
pends on the mass, so that we can estimate that the mini- 
mal mass required for a cluster, in order that the expected 
number of giant arcs in a deep HST observation is ~ 1, is 
~ 2 x 10 14 ft" 1 M G ; 

- the three-dimensional shape of strong lensing clusters does 
not seem to be significantly different from that of the gen- 
eral cluster population. However, strong lensing clusters 
tend to have their major axes oriented along the line of 
sight. This "orientation bias" is larger for clusters which 
produce giant arcs than for clusters which only possess crit- 
ical lines, and it becomes larger by increasing the lensing 
cross section; 

- due to the orientation bias and to the fact that their halos are 
generally well described by prolate triaxial models, strong 
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Fig. 16. The probability distributions of X-ray luminosities of all clusters in the MareNostrum Universe (solid histograms) are 
compared to those of merging clusters (dashed histograms). Each panel refers to different redshift bins, as specified in the title 
above each figure. 
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Fig. 17. As in Fig. 16 but showing the X-ray luminosity distributions of general and merging clusters in the redshift range 
0.25 < z < 0.5. In the left panel we select the clusters with mass 10 14 /z _1 M o < M vil < 5 x 10 14 /j _1 M o . The distributions in the 
right panel refer to clusters with mass 5 x 10 u h~ 1 M Q < M vil < 3 x 10 15 /r'M o (right panel). 



lensing clusters tend to appear rounder when projected on 



the sky. However, zooming over their central regions, we 
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noticed that their projected mass maps are described by 
rather elongated distributions, evidencing the presence of 
substructures projected near cluster cores; 

- the concentrations measured by fitting the density profiles 
of strong lensing clusters stacked in mass and redshift bins 
do not differ significantly from the concentrations of the 
general cluster population. Nevertheless, due to the orien- 
tation bias, the concentrations of the same objects inferred 
from the analysis of the projected density profiles are gener- 
ally larger than in 3D. For clusters with large lensing cross 
sections for giant arcs, the 2D-concentrations can be larger 
by more than a factor of two. These results may provide a 
viable explanation of the large concentrations reported for 
some of the strongest lenses observed so far; 

- the X-ray luminosity-mass relation of strong lensing clus- 
ters is likely to differ from that of the general cluster pop- 
ulation, especially at the lowest masses. We found that at a 
fixed mass, strong lenses with increasingly larger cross sec- 
tions for lensing have higher X-ray luminosities, indicating 
that some process occurring in these objects enhances both 
the X-ray luminosity and the strong lensing cross section; 

- the strong lensing efficiency is certainly correlated with the 
dynamical activity in clusters. We found that clusters with 
large lensing cross sections are characterized by a system- 
atic departure from virial equilibrium. A similar departure 
from virial equilibrium is found for clusters which are in the 
process of accreting substructures, which accidentally tran- 
sit across their cores perpendicularly to the line of sight. 

In conclusion, our results show that strong lensing clus- 
ters are likely to be a very peculiar class of objects, charac- 
terized by several selection biases, which need to be properly 
taken into account in many applications. For example, due to 
the orientation bias, it is very likely that the 3D-masses in- 
ferred from strong lensing models of observed clusters are bi- 
ased high. In arc statistics studies, the statistical modeling of 
the strong lensing cluster population need to include merg- 
ers, triaxiality, and asymmetries in the projected mass distri- 
butions, as also suggested by previous studies. Due to the in- 
trinsic difficulties at modeling analytically all of these effects, 
numerical simulations again seem to be the only viable way to 
describe strong lensing clusters. In this sense, large statistical 
samples of numerically simulated lenses, as that extracted from 
the MareNostrum Universe, are fundamental tools for inter- 
preting the current strong lensing observations. A major effort 
is now necessary to clarify the existing inconsistencies between 
the properties of simulated and observed galaxy clusters, espe- 
cially in the central regions. 
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Appendix A: The most efficient lens in the 

MareNostrum Universe 

We have ranked the clusters in our sample by their lensing cross sec- 
tion, creating a list of the 10 most efficient lenses between z = and 
Z = 2. It turns out that all of them are clusters which are elongated 
along the line of sight and which have significant substructures pro- 
jected near the center. In the upper panels of Fig. |A. 1| we show an 
example, given by the most efficient lens we have identified in the 
simulations. From left to right, shown are the convergence maps at 
three different epochs, namely zi = 0.42, zi = 0.45, and z 3 = 0.48. 
The bottom panels show the cluster at the same epochs as in the upper 
panels, but along a different line of sight, which is perpendicular to the 
previous one. The scale of the figures in the left panels is 633 arcsec. 
In the remaining panels it is 423 arcsec. We can see that: 

1 . in the upper panels the cluster appears rounder and denser than in 
the bottom panels. Clearly, the cluster is elongated along the line 
of sight in the upper plots. In these projections, the lensing cross 
section is larger by more than an order of magnitude compared to 
the projections shown in the bottom panels; 

2. the concentrations measured by fitting the surface density profiles 
of the cluster in the upper panels is significantly larger than that 
obtained by fitting the 3D-density profiles. For example, at z 3 the 
3D-concentration is 4.11, while that inferred from the projected 
mass distribution is 8.29; 

3. going back in time, few substructures approach the cluster center, 
and, at z 3 , they seem to cross the very inner region of the cluster, 
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boosting the lensing cross section significantly. At this epoch the 
lensing cross section is a = 1.92 x 10~ 2 /r 2 Mpc 2 . Between z 3 and 
Z\, it drops by almost a factor of two. At the epoch of maximal 
lensing efficiency the cluster is characterized by rather extreme 
values of the /? and of the r 50 o parameters, which are equal to 
-0.74 and 0.36, respectively; 
4. being the center of the cluster more clumpy, the measured elliptic- 
ity is larger in the central than in the external region. For example, 
at z^ the ellipticities measured within 0.1 x R 2 m and ^ 2 oo are 0.23 
and 0.07, respectively. 

All this shows that the most powerful strong lenses in the universe are 
likely to be a very special class of objects, characterized by several 
peculiarities, which we need to properly take into account in order to 
statistically model them. 

|Oguri & Blandford| (2009) (OB hereafter) used semi-analytic 
methods based on triaxial NFW halos for calculating the probability 
distributions of several properties of the clusters producing the largest 
critical lines in the universe. The size of the critical lines is quantified 
by means of the Einstein radius (see Eq. 18 of OB). Their calcula- 
tions include scatter in concentrations and axial ratios which are cal- 
ibrated with numerical simulations. The length of the critical lines is 
strongly correlated with the lensing cross section (Meneghetti et al. in 
prep.), given that the latter is an area surrounding the caustics, which 
are mapped on the critical lines via the lens equation. Thus, we ex- 
pect that the clusters with the most extended critical lines will have 
the largest cross sections. OB also considered a WMAP1 normalized 
cosmology and adopted several source redshifts for drawing their dis- 
tributions. In particular, they used z s = 1 and z., = 3, which encompass 
the redshift used in our simulations (z., = 2). We attempt now a com- 
parison with their results. 

We begin with the distributions of the cluster orientations and tri- 
axial shapes. They find that in ~ 90% of their Montecarlo realizations, 
the cluster with the largest Einstein radius has 9 < 26 degrees and the 
median of their distribution is between 1 1 and 14 degrees. The cluster 
exhibiting the largest cross section for z., = 2 in the MareNostrum 
Universe has its major axis forming an angle of 23 degrees with the 
line of sight. This cluster has a minor to major axis ratio of 0.21 which 
is also well in agreement with the distributions found by OB, whose 
medians for z s = 1 and z s = 3 are 0.32^ [ and 0.23^^, respectively. 

Considering a WMAP1 normalized cosmology, OB report that 
the typical redshift for the lens with the largest Einstein radius is 
z = 0.28+°-^ for z s = 1 and z = 0.47^ for z s = 3. This is also com- 
patible with the redshift of the cluster discussed above (z = 0.48). The 
distribution of the cluster masses found by OB has a median M vil = 
2.49+°^ x 10 15 /r'M o for z, = 1 andM vir = 1.98!| » x IO^/t'Mq for 
z s = 3. The strongest lens in the MareNostrum Universe has a mass 
of 1.85 x 10 15 /r'M o , which agrees with the results of OB at lcr level. 

As discussed above, the projected ellipticity of the cluster shown 
in the upper right panel of Fig. |A.1| varies between 0.23 in the center 
and 0.07 in the external region. OB find that the cluster with the largest 
Einstein radius has a typical projected ellipticity < 0.3. Given that 
they use simple triaxial mass distributions to model their clusters, the 
projected ellipticity does not change with radius as we find by analyz- 
ing numerically simulated clusters. Nevertheless, our and OB's results 
seem consistent once again. Similarly, there is a very good agreement 
also between the 2D-concentrations: while the most efficient lens dis- 
cussed in this Section has a projected concentration of 8.29, OB find 
that a distribution with median between 8.91^ n 7-®3- 2 ill ls ex ~ 
pected in a WMAP1 cosmology for sources between z s = 1 and z s = 3. 

In conclusions, the most efficient lens in the MareNostrum 
Universe has properties which match the expectations for the clus- 
ter with the largest Einstein radius in the Universe, as derived by 



OB modeling galaxy clusters with triaxial halos, whose structural 
properties are calibrated using numerical simulations different from 
the MareNostrum Universe. This can be considered a valid "sanity 
check" for our results. 
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Fig.A.l. The most efficient lens in the MareNostrum Universe. In the top row of panels we show the convergence maps of 
this cluster at three different epochs, namely at redshifts z\ - 0.42, zi = 0.45, and Z3 = 0.48 starting from the left. The highest 
efficiency for lensing, which corresponds to a cross section for giant arcs of cr — 1.92 x 10~ 2 /T 2 Mpc 2 , is reached at Z3. At this 
epoch, the lensing cross section is approximately two times larger than at z\- The bottom row of panels shows the same cluster 
sequence but along a line of sight perpendicular to that of the upper panels. The side-length of the panels on the left is 633 arcsec, 
while the remaining panels have sizes of 423 arcsec. 



